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For field theories with a topological charge Q, it is often of interest to mea¬ 
sure the topological susceptibility Xt = ((Q^) ~ (Q)^)/^- If we manage 
to perform a Monte Garlo simulation where Q changes frequently, Xt can be 
evaluated directly. However, for local update algorithms and hne lattices, the 
auto-correlation time with respect to Q tends to be extremely long, which 
invalidates the direct approach. Nevertheless, the measurement of Xt is still 
feasible, even when the entire Markov chain is topologically frozen. We test 
a method for this purpose, based on the correlation of the topological charge 
density, as suggested by Aoki, Fukaya, Hashimoto and Onogi. Our studies in 
non-linear a-models and in 2d Abelian gauge theory yield accurate results for 
Xt, which conhrm that the method is applicable. We also obtain promising 
results in 4d SU(2) Yang-Mills theory, which suggest the applicability of this 
method in QGD. 
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1 Motivation 


We are going to address the functional integral formulation of quantum 
physics in Euclidean space with periodic boundary conditions. For a number 
of models of interest, the conhgurations occur in distinct topological sec¬ 
tors, each one characterized by a topological charge Q G Z. Examples are 2d 
Abelian gauge theory and 4d Yang-Mills theories. In these cases also fermions 
may be present, so this class includes the Schwinger model and QCD. Further 
examples are the 0{N) models (non-linear a-models) in iV — 1 dimensions, 
and all 2d CF{N — 1) models. 

We deal with the case where parity symmetry holds, which implies the 
expectation value (Q) = 0. Then the topological susceptibility is given by 

Xt= f d'^x {q{0)q{x)) = , (1.1) 

where q is the topological charge density (Q = f d‘^xq(x)), and Wont is the 
volume. This quantity is often of interest; for instance Xt of quenched QCD 
is relevant for the Witten-Veneziano relation [1]. Clearly, xt can only be 
determined on the non-perturbative level. Hence lattice simulations are the 
appropriate method for this purpose. (Actually the lattice dehnitions of q 
and Q are slightly ambiguous, see e.g. Refs. [2] for comparative studies; we 
will specify later the formulations that we use.) 

Here we consider the Id 0(2) model (quantum rotor), the 2d 0(3) model 
(Heisenberg model), as well as 2d U(l), and 4d SU(2) gauge theories. In our 
Monte Carlo study of non-linear a-models, we apply a cluster algorithm |^, 
which performs non-local update steps. Hence it frequently changes the 
topological sector, so it provides precise results for Xt by direct measurements. 

In most other models of quantum held theory, especially in almost all 
models with fermions or gauge helds, such an efficient algorithm is not known. 
There one resorts to local update algorithms, such as the heatbath algorithm, 
which we used in our gauge theory simulations. In that case, the Markov 
chain tends to get stuck in one topological sector, in particular as one ap¬ 
proaches the continuum limit. That may well happen in lattice QCD with 
light dynamical quarks and a lattice spacing below 0.05 fm |1]. 

In light of these prospects for the near future, indirect methods to measure 
Xt are of interest. Here we test systematically the Aoki-Fukaya-Hashimoto- 
Onogi (AFHO) method [5], which evaluates Xt based on the density cor¬ 
relations {qoqx)\Q\, measured at hxed |Q|. Hence this quantity enables the 
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determination of Xt even from a Markov chain that is entirely conhned to a 
single topological sector. 

An alternative concept with the same motivation is sketched in Ref. |6]. 
For a recent study with a related approach, see Ref. [7]. The procedure of 
Ref. |8] is more general, but it contains yet another option to determine Xt 
from topologically restricted measurements. 

Here we are going to demonstrate in a variety of models that the AFHO 
method works. For suitable settings, it provided Xt values which are cor¬ 
rect to 2 or 3 digits. We will also discuss the practical limitations of this 
approach. 

2 Topological charge density correlation 

The AFHO method was derived in Ref. [5] , inspired by related considerations 
by Brower, Chandrasekhar an, Negele and Wiese [8]. It deals with the long¬ 
distance correlation of the topological charge density at hxed IQI- The 
topological susceptibility Xt can be evaluated from the (approximate) relation 



lim 


( 2 . 1 ) 


where V is the volume in lattice units. The term 



( 2 . 2 ) 


is the kurtosis, which measures the deviation from a Gaussian distribution 
of the topological charges. It tends to be tiny, see e.g. Refs. [9] for quenched 
QCD results, and in the Id 0(2) model it vanishes exactly in the contin¬ 
uum and inhnite volume m- In the current context its contribution can be 
ignored, as we will see in the following. 

Eq. fl2.ip consists of the leading terms of an expansion in 1/(Q^), therefore 
(Q^) = Ext should be large. Since Xt is expected to stabilize in the large 
volume limit, eq. fl2.ip holds up to sub-leading hnite size effects. Moreover, 
its derivation assumes the ratio \Q\/ (Q^) to be small, hence it is favorable to 
apply this method in sectors of small |Q|. 


3 







With these assumptions, eq. fl2.ip shows that the correlation of the topo¬ 
logical charge density in a fixed sector is not expected to vanish over long 
distances. Instead it is expected to attain a plateau, which depends on |Q|: 
it is slightly negative for Q = 0 (obviously, a fluctuation of has to be 
compensated elsewhere), but it rises for increasing Q^. 

The AFHO method has been tested previously in the 2-flavor Schwinger 
model with light chiral fermions m- The numerically measured correlations 
{qo <ix)o suggest that a conclusive evaluation of Xt requires large statistics: on 
a 16 X 16 lattice it requires O(IO^) conhgurations. Variants of this method 
were already applied in 2-flavor QCD, though with a different density [12], 
and recently also in QCD with 2-1-1 flavors, with a reduction to sub-volumes 
[T3] . For a precise test, the non-linear cr-models are perfectly suited, since 
the method can be probed with high statistics, and the results for Xt can be 
compared with reliable direct measurements. In order to probe the potential 
of this approach further, we add investigations in 2d Abelian and 4d non- 
Abelian gauge theories. Synopses of this study were anticipated in proceeding 
contributions [T4] . 


3 Results for the Id 0(2) model 


We start with the Id 0(2) model, or Id XY model, which describes a quantum 
mechanical particle moving freely on the circle 5^, with periodic boundary 
conditions in Euclidean time x. In continuous time, a trajectory can be 
described by an angle (^(x), with (p(0) = Leant)- On the lattice we deal 
with angles (p^,, x = 1... L and (^l+i = (pi|j We introduce the nearest site 
difference 

= {^x +1 - ^x) mod 27r e (-vr, tt] , (3.1) 

i.e. the modulo function is dehned such that is minimized. 

This is one of the simplest models with a topological charge, which is 
given by 


Q 


1 


^-^cont 


dxip'{x) (continuum) , Q 


1 


L 


X=1 


(lattice) , 


(3.2) 


where q{x) = ip'{x)/{27i) is the topological charge density in the continuum, 
and q^ = (27r) is its geometrically dehned counterpart on the lattice. 

^All lattice quantities will be given in lattice units. 
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The continuum action reads S'contl^?] = dxip'{xY. In Appendix 

A we show that relation (12.Ih without the kurtosis term, 

(g(0)g(a:))|Q| =-y + ^ , (3.3) 

is exact in this case, and independent of the separation x. 

In our numerical study, we consider three lattice actions: the standard 
action, the Manton action na and the constraint action [16] . 


5*standard [ 
‘^constraint [ 


L 

COsA^^x) 7 *S*Manton 

x=l 

f 0 A(px < S \/x 
^ +00 otherwise 


^ W A 2 

x=l 

(3.4) 


The parameter /3 (or /dcont) corresponds here to the moment of inertia, and 
in the last case 6 is the constraint angle. 


action 

e 

Xt 

continuum 

2/3cont 

1 

47r2/3cont 

standard 

[in exp{-(3{l-cosip)) 1 

1 <p^exp{-l3{l-cosip)) 

j exp(— ^(1— COS (/?)) COS 9? 

47r^ exp(— jd(l— cos </?)) 

Manton 

b exp(-/3ip2/2) 1 

1 f-T.d.V E exp{-Pip'^ 12) 

^ exp(-/3v52/2)cos(y5 

471-^ exp(-/3!p^/2)) 

constraint 

[ln(5/sin(5))]"^ 

52 

_127r2_ 


Table 1: Closed expressions for the correlation length ^ and the topological 
susceptibility Xt in the Id 0(2) model at inhnite size, in the continuum and 
for three lattice actions. 


In the thermodynamic limit, L —)■ oo, the correlation length ^ = l/{Ei — 
Eq) (he. the inverse energy gap) and Xt are known analytically [IHIIIS], as 
we summarize in Table [To The product ^ Xt is a scaling quantity, i. e. a 
dimensionless term composed of observables. In the continuum it amounts 
to 


iXt 


continuum — 


1 

^ ■ 


(3.5) 


^For the continuum action, Ref. m discusses the spin correlation function, and its 
restriction to a single topological sector. 
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1dO(2), Standard action, p = 2, L = 100 



1dO(2), Manton action, p = 2, L=100 



1dO{2), Constraint action, 8 = 1, L = 100 



Figure 1: The topological charge density correlation in the Id 0(2) model 
over a distance of x lattice spacings, at L = 100. The first two plots refer to 
the standard action and the Manton action at /9 = 2, with ^ = 2.779, (Q^) = 
1.936, and ^ = 4.000, (Q^) = 1.266, respectively. The plot below is obtained 
with the constraint action at 5 = 1, L = 100, ^ = 5.793, (Q^) = 0.844. For 
comparison, we include in all cases horizontal lines for the prediction based 
on eq. fl3.3l) . where we insert the directly measured values of Xt ■ 
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This value is attained for the lattice actions in the limit —)■ oo and <5 —?• 0, 
respectively, which reveals a facet of universality even in one dimension. The 
corresponding scaling behavior is discussed in Refs. [I0 l ll6 p i8] . In particular, 
the Manton action scales excellently, since it is classically perfect. 

Figured] shows examples for numerically measured correlations {cioqx)\Q\^ 
using these actions at L = 100. We see in all cases that the numerical data 
are in excellent agreement with the predicted plateau values. These plateaux 
are accurately visible, so the AFHO method does indeed enable a precise 
numerical determination of Xt- 

To demonstrate this explicitly, we consider the range L = 150 ... 400, and 
IQI = 0, 1, 2, which leads to the results for Xt in Figured For the Manton 
action we obtain precise agreement with the theoretical Xt value in all cases. 
For the standard action we observe small deviations up to a few permille, 
which are suppressed for \Q\ < 1, and for \Q\ = 2 they are reduced as L 
increases. 

These tiny lattice artifacts are revealed due to extremely large statistics: 
for each parameter set, at least 5-10® measurements have been performed 
with a cluster algorithm. This yields very precise results, and illustrates the 
convergence towards the theoretical Xt value for increasing (Q^) (or equiva¬ 
lently L). 

4 Results for the 2d 0(3) model 

We proceed to field theory, and first to the 2d 0(3) model (or Heisenberg 
model), with periodic boundary conditions. In its lattice formulation a clas¬ 
sical spin variable of unit length is attached to each lattice site x, G S"^. 

Regarding the topological charge, we consider sets of three neighboring 
spins. In our case the lattice consists of quadratic plaquettes, and each pla- 
quette is divided into two triangles (the cutting diagonal has an alternating 
orientation between nearest neighbor plaquettes). Each of these triangles 
carries such a set of spins eq, 62 , eq. They are connected on the sphere 5^ by 
the arcs of minimal length to form a spherical triangle. Its oriented area A 
is given by [16] 

X = 1 + Cl • 62 -t- 62 • 6*3 -1- Cl • 63 , ?/ = 63 • (Fa X 63 ) , 

if = arg(x + iy) mod 27r , ^( 61 , 62 , 63 ) = 2ip , (4.1) 

where the modulo function is defined as in eq. fl3.1l) . In each plaquette. 
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Figure 2: The topological susceptibility Xt in the Id 0(2) model for the 
standard and Manton action at /3 = 4, where ^ = 6.8 and 8.0, respectively. 
We show the theoretical value at L = oo, the directly measured value (at 
L = 400), and the values obtained from the AFHO method in the range 
L = 150...400, in the sectors \Q\ = 0, 1, 2. For the standard action, 
there are permille level deviations from the predicted value, in particular for 
IQI = 2, which are suppressed for increasing L. For the Manton action all 
results coincide to an impressive precision, even down to L = 150. 


the normalized sum of these two oriented spherical triangles (the total solid 
angle), = A(ea;, ^a;+i+2)/(47r), is the topological charge density. 

If we sum over all plaquettes, and thus over all triangles, we obtain the 
geometrically defined topological charge 

Q ~ ~ ^ A(ea:, Gx+ 2 ^ ^x+1+2) ^ ^ • 

X X 
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(4.2) 












This definition, which was advocated in Ref. [12] , has the virtue of providing 
integer Q values for all configurations (except for a subset of measure zero), 
just like eq. fl3.2p . It counts how many times (and with which orientation) 
these triangles cover the sphere. 

The standard lattice action reads 


S'. 


standard I ^ I 


X,fl 


(4.3) 


where fi runs from 1 to 2, /i is the unit vector in /i-direction, and f3 > 0. 

Figure [3] shows the topological charge density correlation {qoQx) at /9 = 
1, measured in the sectors |Q| =0, 1, 2, onLxL lattices of size L = 
12 and L = 16. The measurements are carried out parallel to the axes, 
and the spin separation proceeds in steps of two lattice units, due to the 
alternating triangularization in the definition of q^- The horizontal lines are 
the expected plateau values according to eq. fl3.3|) . Again we inserted the 
directly measured values of Xt = they are very precise, thanks to 

the use of a cluster algorithm, which provided a statistics of O(IO^) well 
thermalized and decorrelated measurementsU 

The plots clearly confirm the qualitatively expected picture. We further 
confirm that the condition of a large separation is quite harmless: for the 
separation of four lattice spacings, the plateau value is already well attained. 

For a quantitative analysis, we perform individual fits of the data to a con¬ 
stant in one sector (skipping 2... 6 points at the boundaries); each fit yields 
a value of the topological susceptibility Xt- In addition we consider combined 
fits in two or three sectors. These results are confronted with the directly 
measured values in Figure IH As theory predicts, the lowest |Q|-sectors are 
most reliable. In fact, the evaluation of Xt based on {qoqx)\Q\ is successful to 
an accuracy of a few percent in these cases. However, Figure [3] also shows 
that the application of this method is getting difficult when L increases: then 
the values of {qoqx)\Q\ become tiny, and thus hard to distinguish from zero, 
and from each other. 

This is a case of strong coupling; /3 = 1 leads to a correlation length of 
^ ~ 1.3 (at large L). Hence the volumes that we used can be considered 


^This model is sometimes considered topologically ill, because Xt which is supposed 
to be the scaling term, diverges logarithmically in the continuum limit. In the integral 
representation of eq. (HU, this effect emerges at distance a; = 0; at finite distances, the 
topological charge density correlation is a controlled quantity [16]. Here we determine Xt 
with different methods at fixed so this defect does not affect our study. 
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2d 0(3), Standard action, L=12, (3=1 
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2d 0(3), Standard action, L=16, (3=1 
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Figure 3: The topological charge density correlation in the 2d 0(3) model, 
with the standard lattice action at /3 = 1. We show data measured on 
L X L lattices, L = 12 and 16, in the sectors IQI = 0, 1, 2. They are 
compared to lines for the values according to eq. fl3.3p . with the directly 
measure susceptibility Xt- 


2d 0(3) model, Standard action, p = 1 
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Figure 4: The topological susceptibility in the 2d 0(3) model on 12 x 12 and 
16 X 16 lattices, with the standard action at /5 = 1. We show results obtained 
from hts to the topological charge density correlation in sector |Q|, or com¬ 
bined hts in several sectors. This is compared to the direct measurement as 
(Q ^)/'^ 1 which is very precise. All results are compatible within the errors. 
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large, but the lattice is coarse. In order to probe the AFHO method closer 
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to the continuum limit, we now proceed to a different setting. We move to 
the constraint action, which is defined in analogy to eq. fl3.4p . i.e. the action 
is zero if all angles between neighboring spins are less than 5, and infinite 
otherwise [16]. We set the constraint angle to 5 = 0.55 vr, which corresponds 
to a correlation length of ~ 3.6. Accordingly, we now consider larger 
square lattices, with L = 16 and L = 32. Figure [5] shows the topological 
charge correlations in this case, and Figure [6] displays the fit results. 


2d 0(3), Constraint action, L=16, 5=0.557C 
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-I I I I I I I I I + I I I - 
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Figure 5: Plots analogous to Figure jS] now for the constraint action at 5 = 
0.55 vr, which corresponds to ^ ~ 3.6, and with lattice sizes L = 16 and 


L = 32. 


These plots show again that the condition of a larger separation |x| is not 
a practical problem. However, despite the statistics of O(IO^) measurements, 
the AFHO method does run into trouble in reproducing the directly measured 
Xt values beyond one digit. This is mostly a consequence of the larger volumes 
involved; they suppress the signal, which is relevant to extract Xt in this 
indirect manner. 

Let us finally take a large step to a numerical experiment very close to the 
continuum limit: it is performed with the standard lattice action at /3 = 1.5, 
on square lattices with L = 16 ... 128; at large L, this corresponds to ~ 9.5. 
Figure [7] shows the results for Xt up to L = 84, based on the topological 
charge correlations in the sectors |Q| = 0, 1, 2, and by direct measurement. 
As L increases, the latter converges well at L > 32. The results by the 
AFHO method move towards the directly measured value, and get close to 
it at L = 40. Here the range L ^ 40... 60 is optimal for its application. 
As we increase L further, we face again the problem that the tiny signal. 
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2d 0(3) model, Constraint action, 5 = O.SSti 
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Figure 6: The topological susceptibility in the 2d 0(3) model on 16 x 16 and 
32 X 32 lattices, with the constraint action at 5 = 0.55 tt, where ^ ~ 3.6. 
We display results from hts to the topological charge density correlation in 
sector IQ I, or combined fits in several sectors. Comparison to the direct 
measurement shows that the AFHO method is less successful than in the 
examples of Figure 01 since larger volumes are involved. 

which matters for Xt, gets lost in the statistical noise. For |Q| = 1, 2 this 
happens already at L > 64; only at Q = 0 the method still leads to useful 
results up to L = 84. (For completeness we add that at L = 128 we obtained 
Xt = 0.0019(8), which is still compatible with the directly measured value 
0.002292(9), but it has an error of 42 %). 

5 Results for 2d Abelian gauge theory 

We proceed to 2d U(l) gauge theory with the plaquette action, i.e. Wilson’s 
standard lattice formulation [20]. We simulate it with the heatbath algorithm 
on periodic L x L lattices, with a variety of L and {3 values, which keep 
(Q^) in the range of 0.7 to 10.4. In each case, the statistics involves 10^ 
conhgurations. For the topological charge density we also applied the 
straight plaquette regularization of the held strength tensor in terms of non¬ 
compact link variables 
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2d 0(3), Standard action, (3 = 1.5 
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Figure 7: The topological susceptibility Xt for the standard action at /3 = 
1.5, on L X L lattices with L = 16... 84 (with ^ ~ 9.5 at large L). The 
directly measured values stabilize for L > 32 (its errors are too small to be 
visible in this plot), and the AFHO results approximate it well in the regime 
L = 40 ... 84. For smaller L, this method suffers from significant finite size 
effects, and for larger L the signal for the determination of Xt is too small 
for a good numerical resolution. 


still with the modulo function as defined in eq. iu. As in Section 4, its 
correlation was measured parallel to one of the axes. 

Thanks to a generous separation of the measurements, the direct evalua¬ 
tion of Xt is very precise, although the updates are local. 

Figure E] gives three examples which illustrate that also here the data for 
the topological charge density match the predicted plateaux well, so that the 
determination of Xt by the AFHO method is possible. In addition we see 
again the difficulty setting in as the volume increases. 

Figure [9] provides an overview over the results in the range L = 6... 20 
and /3 = 1... 5. In the extreme cases, the plaquette values amount to (P) = 
0.4464 (/3 = 1) and 0.8934 (/? = 5); for the rest we refer to the caption of 
Figure [HI 

In this case, we apply the AFHO method as a combined fit to the data in 
the topological sectors with Q = 0 and [Ql = 1- This might be the optimal 
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2d U(l), L=16, p=3 
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Figure 8: The topological charge density correlation in 2d U(l) gauge theory 
on L X L lattices, over a distance of x lattice spacings. We show the hrst 
three plateau values (horizontal lines, based on direct measurements of Xt), 
and the corresponding data for (3 = 2, L = 12, with a plaquette value of 
(P) = 0.6978 and (Q^) = 2.79; (3 = 3, L = 16, (P) = 0.8100, (Q^) = 2.83, 
and /3 = 4, P = 20, (P) = 0.8635, (Q^) = 3 . 02 . 
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application, which includes a large number of data points, but avoids the 
less reliable topological sectors. In all the cases shown in Figure [HI these 
AFHO results for xt agree with the directly measured values, within errors 
on the percent level, e.g. at L = 16 we obtain Xt = 0.0196(6) {/3 = 2); 
0.0110(3) {/3 = 3); 0.0075(2) {/3 = 4). 
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Figure 9: The topological susceptibility Xt for 2d U(l) gauge theory, on 
L X L lattices. The horizontal lines are the directly measured values (errors 
and differences for various L are invisible on this scale). The symbols are the 
AFHO values obtained by a combined £t in the sectors |Q| < 1. 

In contrast to the previous sections, we are now dealing with a local up¬ 
date algorithm, which is the situation that motivates this project. As an 
illustrative example, we add a measurement of the integrated autocorrela¬ 
tion time Tint (for the dehnition, see e.g. Ref. [2T]) with respect to Q. It is 
expressed in the number of sweeps (updates of each link variable), in a fixed 
volume with L = 16. Figure fTOl shows that Tint increases rapidly as we ap¬ 
proach the continuum limit. This conhrms that, for the heatbath algorithm, 
the problem of “topological freezing” becomes severe indeed, once we attain 
plaquette values of (P) > 0.9. (The results in Figures E] and [9] were ob¬ 
tained by separating the measurements by a number of sweeps, which clearly 
exceeds Tint). 
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2d U(l), L=16 



P 

Figure 10; The integrated autocorrelation time Tint in 2d U(l) gauge theory, 
with respect to the topological charge Q. We consider = 1... 5 on a 16 x 16 
lattice, and measure Tint based on 10^ sweeps, after thermalization. As (3 is 
getting large, we observe an exponential increase of Tint (measured in the 
number of sweeps). 

6 Results for 4d SU(2) Yang-Mills theory 

Finally we study 4d SU(2) Yang-Mills theory, which is in many respects 
close to QCD, but computationally cheaper. Therefore, tests in SU(2) Yang- 
Mills theory are ideal for hints about the practical applicability of the AFHO 
method to QCD. As usual, the lattice gauge field is represented in terms of 
compact link variables G SU(2) [20] . 

Following Ref. [22], we determine the topological charge density and 
the topological charge Q = Yhx improved field theoretic definition, 

= Y YY’lt'l • (6-1) 

□ = 1,2,3 

where denotes the dimensionless lattice held strength tensor, clover 

averaged over square-shaped Wilson loops of size □ x □, and Ci = 1.5, C 2 = 

—0.6, C3 = 0.1. 

We apply eq. (I6.1|l . after performing a number Ncooi of cooling sweeps 
with the intention to suppress UV huctuations in the gauge conhgurations. 
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while preserving the topological structure. A cooling sweep amounts to a 
local minimization of the action, i.e. a successive minimization with respect 
to each gauge link. For this minimization we also use an improved version of 
the lattice Yang-Mills action, 

= 4 Y Y E ^Tr(l - , (6,2) 

X fiu □=1,2,3 

where (3 = and is a clover averaged loop of size □ x □. 

Choosing an appropriate sweep number Ycooi is a subtle and somewhat am¬ 
biguous task, which will be discussed below. 

The lattice action used for the generation of gauge configurations is the 
standard plaquette action, which is obtained from eq. fl6.2p by setting ci = 1, 
C2 = C3 = 0. As in Section 5, the simulations were performed with a heatbath 
algorithm [20], now at /3 = 2.5. This corresponds to the lattice spacing 
a ~ 0.073 fm, when the scale is set by identifying the Sommer parameter rg 
with 0.46 fm [23]. That value of a is in the range of lattice spacings 0.05 fm < 
a < 0.15 fm typically used in contemporary QCD simulations. We generated 
about 4000 configurations in each of three volumes, V = 14“^, 16^ , 18“^. This 
is also a typical statistics in QCD simulations. 

Assigning a topological charge Q to each configuration leads to the statis¬ 
tics given in Tabled for the sectors |Q| < 4. We proceed by computing the 
correlation function of the topological charge density (goQx}iQi in all these 
sectors Q and volumes V. 

The normalization factor on the right-hand-side of eq. fl2.ip is given by 
the inverse volume. As we have observed in the previous sections, the corre¬ 
spondingly suppressed signal in a large volume is often the bottleneck in the 
application of the AFHO method. In order to compensate this suppression, 
which is worrisome in a 4d volume, we now determine (qoqx) by measur¬ 
ing all-to-all correlations in each configuration, thus taking advantage of the 
discrete translational and rotational invariance. 

In a second step we fit to these lattice results the right-hand-side of eq. 
fl3.3p . i.e. eq. fl2.ip with C4 = 0, with respect to Xt, at sufficiently large sepa¬ 
rations X, where {qoqx)\Q\ exhibits a plateau. Figure ITT] illustrates the deter¬ 
mination of Xt after Ycooi = 10 cooling sweeps, in the three lattice volumes 
under consideration. Clearly the correlation function {qoqx)\Q\ is different 
for each topological sector |Q|. These differences are more pronounced for 
smaller volumes V and larger topological charges Q, which shows that the 
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V 

Q = 0 

IQI = 1 

IQI = 2 

IQI =3 

IQI=4 

144 

1023 

1591 

893 

350 

103 

16^ 

722 

1371 

942 

574 

248 

18^ 

622 

1079 

898 

616 

402 


Table 2: Number of configurations for three volumes V, in each topolog¬ 
ical sector 0 < IQI < 4. The topological charge has been assigned after 
performing A^^cooi = 10 cooling sweeps. 


splitting according to eqs. fl2.ip and fl3.3p can indeed be resolved from our 
data. Hence the statistics, drastically amplihed by the all-to-all correlations, 
is indeed sufficient to reveal the relevant signal. 

In particular, for V = 14^ the maximally available on-axis separations 
|x| = 6, 7, 8 are at the border-line which allows us to observe plateaux of 
(go<?x)|Q|- For V = plateaux are visible in the range |x| = 7,8,9 and in 
V = 18^ even hve points, 7 < |x| < 11, are consistent with a plateau. 

Since the signal, i. e. the differences between the plateau values, increases 
for smaller volumes, a promising strategy might be to use anisotropic vol¬ 
umes. For example the volumes of a 14^ and a 12^ x 24 lattice are similar 
(he. both should exhibit a similar signal quality), but the latter allows to 
study larger separations (on-axis up to |x| = 12, though not with the entire 
statistics of all-to-all correlations). 

In Figure [12] we show determinations of Xt and compare different numbers 
of cooling sweeps, Ncooi = 5 , 10 , 20, in the volumes V = 14*^ and 18^^. For a 
small number, such as Ncooi = 5, the correlation function {qo(lx)\Q\ is rather 
noisy. This is a consequence of strong UV fluctuations, which are manifest in 
the topological charge density and which are not hltered out sufficiently at 
small A^^cooi- For a large number of cooling sweeps, like A^cooi = 20, statistical 
errors are signihcantly smaller, but the correlation function {qoqx)\Q\ exhibits 
plateaux only at larger separations |x|. 

This effect becomes plausible when considering the structure of the states 
contributing to {qoqx)\Q\- This correlation function is the Fourier transform 
of an analogous correlation function, summed over all topological sectors 
at a hnite vacuum angle 9 (for a detailed discussion, see Ref. [H]). The 
plateau values arise due to the non-vanishing vacuum expectation value (q^) 
at 6* 7 ^ 0. Deviations from these plateaux are predominantly caused by 
low-lying excitations, which correspond to glueballs in Yang-Mills theory. 
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4dSU(2),U14, N,„„|=10 



|x| 

4dSU(2),L=16, N,„„|=10 



|x| 

4dSU(2),L=18, N,„„|=10 



|x| 


Figure 11: The correlation {qQqx)\Q\ as a function of the on-axis separation |x|, 
after N^oox = 10 cooling sweeps, for the lattice volumes V = 14*^, 16^, 18"^. 
Fits of the right-hand-side of eq. fl3.3l) with respect to Xt are indicated by 
the horizontal solid lines. 
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Figure 12: The correlation {qQqx)\Q\ as a function of the on-axis separation 
|a:| for different numbers of cooling sweeps, Ncoo\ = 5, 10, 20, and lattice 
volumes V = 14^ and 18"^. 


Due to the glueball size, the overlap with ga,]!!) (where ID) is the vacuum 
state) increases when using extensive cooling (then q^ is an extended operator 
resembling a low-lying glueball), compared to little or no cooling (then q^ is 
a highly local operator). Consequently, cooling enhances the contribution of 
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excitations to the correlation fnnction {qoqx)\Q\ and, hence, canses stronger 
deviations from the plateanx. 

In practice one shonld decide for an optimal compromise, i.e. an inter¬ 
mediate nnmber of cooling sweeps. Snch a compromise can be read off for 
instance from plots showing the dependence of the topological snsceptibility 
Xt (obtained by the AFHO method at hxed |Q|), or the correlation fnnction 
{qoqx)\Q at specihc separation |a;| ~ L/2 on Ncooh Snch plots are shown in 
Fignre[T3]for L = 18 and hxed topological charge |Q| = 0,1, 2. As expected, 
the statistical errors are qnite large for small A'^cooi ~ 5. In an intermediate 
region, 5 < A^^cooi S; 15, there are stable plateanx of both {qoqx)\Q\ and Xt with 
comparably small statistical errors. The rather long plateanx indicate that 
cooling is a nnmerically stable procednre, not destroying topological excita¬ 
tions nor introdncing any nnwanted non-locality effects. At large A^"cooi ~ 15 
there is a slight trend towards lower xt, in particnlar for Q = 0, which conld 
be a hrst sign of contamination by excited states. An optimal choice for A'^cooi 
is somewhere inside the plateanx region, snch as A^^cooi = 8 or A^cooi = 10, as 
we have nsed in the examples in Fignre [TT] and Table [31 

Nnmerical results for the larger volumes, V = 16"^ and 18^, and moderate 
cooling, A^^cooi = 8 or 10, where a reasonably accurate determination of Xt 
seems possible, are summarized in Table [3l Fits have been restricted to sec¬ 
tors Q fulhlling \Q\/{xty) < 0.5, which involves the small expansion parame¬ 
ter in the derivation of eq. fl2.ip [MT7] . The resulting values for the topological 
susceptibility agree within the errors with a previous straight determination 
(without topology hxing), which arrived at Xt = 7.0(9) x 10“® |22] . 


V 

^coo\ 

Q = 0 

IQI = i 

IQI = 2 

IQI = 3 

combined 

16^ 

8 

6.1(9) 

5.5(9) 

7.1(8) 


6.3(6) 

6.3(5) 

10 

7.4(10) 

5.4(7) 

6.5(9) 


18^ 

8 

7.1(14) 

5.8(10) 

7.4(12) 

9.2(19) 

7.3(6) 

7.0(5) 

10 

6.2(10) 

5.9(10) 

6.6(9) 

8.7(11) 


Table 3: Results for the topological susceptibility Xt x 10^ extracted from 
hts to a single sector [Ql, or combined hts to several sectors (0 < [Ql < 2 for 
V = lO'^, 0 < IQ I < 3 for R = 18"^). A corresponding study without topology 
hxing [22] obtained Xt x 10^ = 7.0(9). 

We also investigated the magnitude of ordinary hnite volume effects, (not 
related to topology hxing). We did this by computing and comparing Xt at 
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4d SU(2), L=18 



4d SU(2), L=18 



Figure 13: The correlation {qoqx)\Q\ at |x| = L/2 = 9, and the topological 
susceptibility Xt, as a function of the number of cooling sweeps A^"cooi in the 
volume V = 18^. 


unhxed topology using Xt = {Q‘^)/V- For the three volumes used throughout 
this section, and N^ooi = 10, we obtain for XtiY) 

Xt(14") = 7.02(1) X 10-^ xt(16^) = 6.97(2) x 10-^ Xt(18^) = 7.01(3) x 10"^ 

These values agree within the statistical errors. Moreover, both the differ¬ 
ences as well as the statistical errors are signihcantly smaller (by more than a 
factor of 10) than the uncertainties associated with our xt determinations us¬ 
ing the AFHO method, as listed in Table |3l Therefore, in our study ordinary 
hnite volume effects can safely be neglected. 

The similarity of 4d SU(2) Yang-Mills theory and QCD suggests that an 
application of the AFHO method to QCD will also allow for a determination 
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of the topological susceptibility up to about 10 %, based on 0(1000) configu¬ 
rations for each of the topological sectors with small IQI- We plan to explore 
this scenario in the near future. 


7 Conclusions 

We have investigated the AFHO method [5] for the evaluation of the topo¬ 
logical susceptibility Xt based on the correlation of the topological charge 
density. Amazingly, this method allows — in principle — for a measurement 
of Xt even within a hxed topological sector, and low |Q| are most promising. 

We have seen that the method as such works; in some cases it provided 
results, which are correct to 2 or 3 digits. Hence situations do exist, where 
the approximations in the derivation of formula fl3.3p (in particular including 
order 0(1/1/^) incompletely, and neglecting all terms of 0{1/V^)) can be 
justihed, as we could conhrm on the non-perturbative level. 

This approximation neglects the kurtosis term in relation fl2.ip . and all 
higher terms. They are suppressed by the inverse volume. Moreover, it 
is a generic property that topological charges are approximately Gauss dis¬ 
tributed, such that |c4| is small. As a peculiar case, it vanishes in the contin¬ 
uum Id 0(2) model, and in that case there are no higher corrections either 
(cf. Appendix A). 

Regarding the limitations of the applicability range, we note that — for 
most models studied here — the theoretical condition of measuring the cor¬ 
relation “at large separation” turned out not to be worrisome in practice; 
only in Section 6 this issue has some relevance. However, the AFHO method 
generally runs into trouble when the volume V increases. Then the signal 
in (goQx}iQi is suppressed by a factor 1/H, and the separation between the 
predicted plateaux by 1/H^. If we want the correlation length to be clearly 
larger than the lattice spacing (so that lattice artifacts are under control), 
we need a sizable lattice to keep the hnite size effects under control as well. 
With these conflicting requirements, even in our 2d test models, and despite 
a statistics of 10^ measurements, this method led to results with rather large 
statistical errors. 

Thus we observe that the AFHO method, which is based on topologi¬ 
cally restricted measurements, is plagued by unusually persistent hnite size 
effects. In usual settings, the latter are exponentially suppressed, i.e. they 
are oc exp(—const. L/^), and quite small if L/^ > 4. However, topologically 
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restricted numerical measurements are very sensitive to finite size effects: we 
recall that relation fl2.ip is a truncated polynomial expansion in 1/ld, but 
we can also refer directly to the large ratios L/^, which were used all over 
this study: e.g. in the upper plot of Figure |2] it amounts to L/^ > 22, which 
usually makes hnite size effects negligible, but here they are signihcant. This 
property calls for a larger size L. In turn, that causes problems in extracting 
the subtle effect, which is relevant for the indirect determination of Xt- 
That issue might be the bottleneck for the prospects to apply the AFHO 
method in four dimensions, and it has inspired the approach of restriction to 
sub-volumes [13] . However, our study in Section 6 shows that the suppression 
of the wanted signal by the inverse volume can be successfully compensated 
if we enhance the statistics by means of all-to-all correlation measurements. 
We have seen in 4d SU(2) Yang-Mills theory that this procedure does pro¬ 
vide sufficient precision for {qoQx), even for a moderate number of 0(1000) 
conhgurations in one topological sector. Due to this observation, the AFHO 
method in its original form appears quite promising in QCD. We are going 
to test if it enables also there the evaluation of Xt in an unconventional way, 
to an accuracy of about 10 %. 

Acknowledgements Lilian Prado has worked on this project at an early 
stage. We also thank Christopher Czaban and Philippe de Forcrand for 
interesting discussions. This work was supported in part by the Mexican 
Consejo Nacional de Ciencia y Tecnologm (CONACYT) through project 
155905/10, by DGAPA-UNAM, grant IN107915, and by the Helmholtz In¬ 
ternational Center for FAIR within the framework of the LOEWE program 
launched by the State of Hesse. A.D. and M.W. acknowledge support by 
the Emmy Noether Programme of the DFG (German Research Foundation), 
grant WA 3000/1-1, and G.P.H. acknowledges support through the project 
Redes Temdticas de Colaboracion Academica 2013, UGOL-GA-56. Galcula- 
tions on the LOEWE-GSG and on the FUGHS-GSG high-performance com¬ 
puter of the Frankfurt University were conducted for this project. We would 
like to thank HPG-Hessen, funded by the State Ministry of Higher Education, 
Research and the Arts, for programming advice. 


24 



A Topological charge density correlation in 
the continuum Id 0(2) model 


The key quantity of this work is the correlation function of the topological 
charge density. In this appendix we compute it analytically for the Id 0(2) 
model, formulated in continuous Euclidean time x, with periodicity length 
L. For this purpose, it is useful to include a 6*-term in the action, 

Q ^Tcont 

= dx^\xr-i9Q[^]. (A.l) 

In the canonical formulation of quantum mechanics, the corresponding Hamil¬ 
ton operator, its energy eigenfunctions and eigenvalues read m 


H 


1 

2/3cont V 271) 


{ip\n) 



1 ( 

2/3eontV 277 ) ’ 


(A.2) 


where p = —i^ and n E h. The operator for the topological charge density 
is given by 


= —[H,0] = —^— 
27r^ 27r/3cont 


dp 271) 


(A.3) 


This operator is anti-Hermitian (due to the Euclidean time derivative of the 
Hermitian operator p), with the matrix elements 


(m|g|n) = 


(27r)2/3cont J-n 
Hence the expectation value 

{q{x)) = ^ 


dp e 


—imip 




dp 


271J 


(2vr)2/3, 


cont 


Z{e) {277f^oni 


— 27in) 


) — EnLcont 


, with Z{9) = 


3 £/n,Tcont 


n£’L 




is imaginary in general (of course it vanishes at 0 = 0). 

Therefore the corresponding correlation function is in general negative^ 


(g(O)g(x)) = {m\q\n) {n\0m) e 


EnX Em (Tcont 3?) 


1 1 


Z{9) {277Yf5f 


cont 


^(0 - 27m) 


2 g ^nTcont 


(A.4) 


n^lj 


field theoretic models, the correlation of the topological charge density over large 
distance is known to be negative as well [24]. 
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It is remarkable that this correlation is independent of x, if x/L^ont ^ ^ (this 
condition allows us to insert a unit factor between the end-points). 

The vacuum angle 6 enables also the computation of the topologically 
restricted partition function, 


^Q = 


27r 


dd z{e) _ 


,-QV(4a) 


a = 




cont 


and correlation function, 

(?(0)g(a:))Q = 


27rZ 


Q J — 7 \ 

1 


dO Z{e){q{fd)q{x)) 


PlontZq y_c 

1 




- 1 + 


doe^e 




2 -ae'^-ieq 


2a 


Finally we insert a and Xt = tt-f^cont/2 (cf. Table [T]) to arrive at 


(A.5) 


{q{0)q{x))Q = -j^ + .5— . (A.6) 

.^cont -^cont 

Also the topologically restricted correlation function is constant in x, which 
explains that the data in Section 3 attain the plateau values immediately. 

Moreover, we see that eq. (13.31) is exact in this specihc case, which is 
consistent with the fact that the kurtosis vanishes [10]. Therefore, in our 
numerical study presented in Section 3, the actual issues are lattice artifacts 
and the visibility of the predicted plateau values in numerical simulation 
data; both are generally relevant questions. 
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